CRITICAL BEHAVIOR OF A FOUR-POINT 

SCHWINGER-DYSON EQUATION 



We study the Schwinger-Dyson equation associated with a chirahty- 
changing fermion 4-point function in a strongly-coupled U{1) gauge theory. 
After making appropriate simplifications, we solve the equation numerically 
via a relaxation method. Our analysis provides an estimate of the critical cou- 
pling and it gives some indication as to the general momentum dependence 
of the 4-point function. 
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1 Introduction 



In this paper we wish to study the dynamical generation of momentum dependent 
fermion 4-point functions within strongly interacting gauge theories. These will be 
purely nonperturbative quantities associated with breakdown of chiral symmetries, 
and would for example imply the existence of the corresponding 4-fermion conden- 
sates. Our analysis is based on the Schwinger-Dyson (SD) formalism. The resulting 
equations are analytically intractable and thus we seek numerical solutions. 

We have two main goals in this study. One is to estimate the "critical coupling" 
necessary for the formation of the 4-fermion condensates. We would like to study 
this in the context of a situation where 4-fermion condensates develop on scales 
higher than the scale of mass formation. This would make it consistent to treat the 
4-fermion condensate problem in isolation, independently of the mass generation 
problem. (This is not the case for QCD, but even in that case it is useful to at 
least consider the existence of the effects we study here). In extensions of the 
standard model the hierarchical symmetry breaking pattern we envisage can be a 
natural consequence of the fact that a 4-fermion condensate may break fewer gauged 
symmetries than the 2-fermion condensate. 

For example, in [1] it is shown how the top mass may be associated with a cer- 
tain 4-fermion condensation occurring well above the electroweak symmetry break- 
ing scale. The point is that a 4-fermion condensate, composed of fermions with 
standard quantum numbers, can signal a dynamical breaking of isospin symmetry 
while not breaking SU{2) x U{1) gauge symmetry. A pattern of symmetry breaking 
is described [1] in which the isospin breaking feeds strongly into the top mass while 
not significantly infiuencing the W and Z masses. 

Our other goal is to make a first attempt to extract the momentum dependence 
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of the dynamical fermioii 4-point function. For example if such an object is to 
play a role in generating a fermion mass, then two lines of the 4-point function 
may be closed off into a loop. The details of the momentum dependence, such as 
knowing the relative size of the 4-point function when different pairs of momenta 
are large, is then of interest. As a first guess, it is standard to appeal to some kind 
of factorization or vacuum saturation and treat the 4-point function as a product of 
two 2-point functions. We proceed beyond that type of picture here, although we 
are forced into making our own truncations of the full treatment. 

Truncations are also made in the study of the dynamical mass function, but 
there is growing confidence that most of the essential physics of chiral symmetry 
breaking is contained in these treatments. For example, the kernel of the SD equa- 
tion is usually dealt with in some form of ladder approximation. But in the case 
of a constant coupling, which is the case we study here, it is found that the main 
conclusions regarding the behavior of the mass function are quite insensitive to the 
precise form of the kernel [2] . 

In the present study we use a one gauge-boson exchange approximation. Since 
the gauge boson can attach to any pair of the four legs, the SD equation sums up 
a much more complicated set of diagrams than the set of ladder graphs. It is clear 
from the study of the mass function that the poorly understood infrared effects (e.g. 
confinement) are not important in extracting the intermediate and high momentum 
behavior of the mass function. We shall adopt a similar philosophy here, and drop 
various nonlinearities which arise only in the infrared. 

In the next section we consider the SD equation for the fermion 4-point function. 
In section 3 we present an explicit derivation of the equation within the framework 
of certain approximations, which are further discussed in section 4. In section 5 we 
discuss our numerical treatment, and in section 6 we present our results. 
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2 Schwinger-Dyson equation 

Our model consists of a strongly interacting U(l) gauge theory with one Dirac 
fermion. We start by considering 4-fermion operators of the form '^gT abi^bi^J^'cd^d^ 
where F and F' are 4x4 matrices and a, h, c, d are spinor indices. In this work, we 
are particularly interested in operators that are purely chirality changing^ i.e. of 
the form ipLTilijiipLr'ipji + h.c, where i/jl = ^^-^'0 i^R = ^^^^i^- The two 
independent operators which have this property and respect parity are^ 

ija^'^iljija^^xp = il)L(T^''il;RiljLa^''ijR + h.c.. (1) 

The vacuum expectation values of these two operators receive exclusively non- 
perturbative contributions, and their non-zero values signal a dynamical break-down 
of a chiral U{1)a- However the discrete symmetry e*^"/^^'^'''^ where n is an integer 
can remain unbroken and prevent fermion masses. We imagine that this discrete 
symmetry is broken on a smaller energy scale in which case the resulting fermion 
masses can be neglected in our analysis. 

We develop the appropriate SD formalism in analogy with the standard anal- 
ysis of the 2-point function. The SD equation for the 2-point function expresses 
this function in terms of a 3-point function, as shown diagrammatically in Fig. 1. 
This equation is part of an infinite SD hierarchy, which must be truncated in some 
manner. The loop expansion in the elTective potential formalism [4] leads to such 
an approximation, and at two-loops it yields the integral equation given in Fig. 2. 

The diagrammatic form of the SD equation for the 4-point function is shown 

in Fig. 3. It involves a 5-point function, and it is again a part of the infinite SD- 
^For a discussion on chirality-conserving 4-fermion operators, see [3]. 

^We choose to work in Euclidean space, so there is no difference between upper and lower 
Lorentz indices. 
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equation hierarchy. In the corresponding effective potential we have terms as shown 
in Fig. 4a, b. The smaher black dots depict the fermion self-energies, and the larger 
dots the 4-point function. If the fermions have non-zero self-energies, terms of the 
form shown in Fig. 4a, with and without gluons attached, would be important. 
The 4-fermion condensate would turn out to be proportional to the square of the 
2-fermion condensate, if other diagrams are neglected. 

However, if the fermion masses can be neglected, as we assume in this work, then 
the diagram in Fig. 4b is essential for a nonvanishing 4-point function. Minimization 
of the effective potential with respect to the 4-point function leads us to an equation 
of the form shown in Fig. 5. It is clear that the 5-point function of Fig. 3 has been 
approximated by a 4-point function, with a gauge boson attached to one of the 
fermion legs with a bare coupling, as shown in Fig. 4c, and in complete analogy 
with the study of the 2-point function. Note that we have symmetrized the set of 
diagrams such that a gluon is attached to all possible pairs of fermion legs. We then 
have to insert a factor of 1/2 to recover the correct SD equation. 

The resulting SD equation is linear. Non-linearities would come from vacuum 
terms like the one shown in Fig. 6a, which contributes to the SD equation a term 
displayed in Fig. 6b. Note that the number of black dots in Fig. 6a must be even. 
These non-linear terms are not easily resummable, in contrast to the case of the 2- 
point function or the Nambu-Jona-Lasinio model. We do not consider these terms in 
this work since we expect them to affect the 4-point function mostly in the infra-red 
regime, in analogy with the 2-point function. 

We now derive the explicit analytic form of the equation shown diagrammati- 
cally in Fig. 5. 
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3 Feynman diagrams 

We are considering a 4-point function associated with the Green's function 
0\T{ijja{x)ipb{y)ipc{z)ipa{0)}\oY The 4-point function receiving exclusively non- 
perturbative contributions is, in momentum space 

where pi, ...,p4 are the external momenta corresponding to the fermions with spinor 
indices a, ...,d respectively. We wish to develop the SD equations for the two scalar 
functions Os+p and Or- 

When these functions appear inside the integrals in the SD equation, then their 
arguments depend on the loop momentum k flowing in the gauge boson propagator. 
For example, we have the k dependent variables s\f, and u\ where the superscript 
i is equal to A, ...,F, corresponding to each of the six diagrams on the right-hand 
side of the equation in Fig. 5. They are given as follows:^ 

s = -iPi+P2y t = -{pi-pif 

= s = -{pi - Pa - 

= s = -{pi- pi-\-kf 

= -(pi + p2 - kf t'^ = t 

= -iPi+P2 + kf t^ = t 
= = 
= = t^ 

The 4-point functions with arguments as they would appear inside the six diagrams 
are denoted by C^+p and C^, where i = A, F. 



u = - 


-{pl 


= 


-ip 






u"" = 








u'^ = 








u^ = 




= 


u 


u^ = 


u, 



^Wc define those variables with a minus sign because we work in Euclidean space, and we want 
to work with positive quantities. 
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We then define the following functional operators F*: 

r^lK] = —J d'k- ^ 



Stt^ J {pz + kf{pi-kf 
T^\K\ = A / d^k. — 

r^[K] = A / d^k. — 

r^[K] = A / TT^ 7T^ 

where K is a function of the loop and external momenta, with a possibly non- 
trivial spinor structure, and a is the momentum-independent coupling. We work 
in the Landau gauge which is popular in studies of SD equations; the gauge boson 
propagator reads ^ = ^ (^6^" - 

With regards to the gauge dependence of our results we note that similar gauge 
dependence occurs in the two-point SD equation, and that the choice of Landau 
gauge in that case gives a critical coupling resembling that found in more complete 
treatments [5], [6]. Since one of our main objectives is to compare the critical 
coupling for the two-point and four-point functions, it seems natural to use the 
same approximation (ladder) and the same gauge (Landau). We also note that in 
our case of massless fermions the one-loop fermion self-energy vanishes in Landau 
gauge, thus making the use of the bare photon-fermion vertex consistent with the 
Ward-Takahashi identity. To do better and to approach our problem in a truly 
gauge invariant way, we would have to carry over the sort of treatment presented in 
[5]; this is a task far beyond our present means. 

To significantly reduce the complexity of our problem we make the following 
truncation. Inside the loop integrals we drop certain terms which are suppressed in 



the large loop-momentum limit. In particular, we drop terms proportional to , 
i = 1...4, coming from the fermion propagators, and we replace k'^k" terms, where k 
is the loop momentum, by ^8^^ . An example of the effect of this type of truncation 
in a SD equation can be found in the literature. The authors of Ref.[7] used a 
3-point SD equation for a calculation of the pion decay constant /jr. They found 
that results obtained by neglecting terms in the pseudoscalar vertex proportional to 
external momenta were surprisingly very similar to their full results. The truncation 
in that case turns out to correspond to the popular Pagels-Stokar approximation 
[8] for the calculation of Our case is of course quite different, and we will test 
further the effect of such a truncation below. 

In the following, by "Scalar insertion" we mean placing the function Os ■ I ® I 
inside the integral on the right-hand side of the SD equation, by "Pseudoscalar 

The spinor structure can be simplified as follows: 
I. Diagrams A and B 

a) Scalar insertion: 



insertion" placing Op • 7^ (8) 7^, and by "Tensor insertion" placing Ot ■ cr^^ ® cr^^ 



o 



b) Pseudoscalar insertion: 



o 



= 3r^'^[C^'^]75®7 



c) Tensor insertion: 



YA,b 



o. 



a 



TP 



a- 



II. Diagrams C and D 
a) Scalar insertion: 



^C,D 



b) Pseudoscalar insertion: 
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c) Tensor insertion: 



where we used the identity 
III. Diagrams E and F 



7^) 



al Scalar insertion: 



YE,f 



o 



b) Pseudoscalar insertion: 



cl Tensor insertion: 



o. 



We may now equate terms proportional to /(g/, 7^(g7^ and a^'' ®a^'' on either 
side of SD equation, in order to obtain the following system of coupled integral 
equations: 

Os+P = 3{r^ + r^)[Os+p] + 

6(r^ + r^-r^-r^)[0r] 
Ot = (-(r^ + r^) + 2(r^ + r^ + r^ + r^)) [Or] 

+l(fc + r^_r^_r^) ps+p] , (4) 

where we have defined = r[0'],i = A,...,F. 
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4 Further simplification 

This is still a formidable system of integral equations when it is realized that the 
quantities Og+p and Ot are functions of six variables, s, t, u, pi, pi, and p^. The 
most interesting dependence is on the three variables s, t, and u, and we would be 
happy if we were able to solve for the function 



where = s + t + u. But the SD equation does not permit this, since for loop 
momenta much larger than the external momenta we must know the value of O 
with two of the pi's much larger than the other two. 

In order to proceed we shall make the following replacement inside the integrals: 



A prime indicates a possible dependence on the loop momentum k. For large loop 
momenta, we find it reasonable to assume that the value of O lies between its value 
when the large momentum flows through all of its legs and its value when the large 
momentum is removed. We have chosen the geometric mean. 

This prescription yields equations which can now be solved for the function 
0(s,t,u). The final equations are 



dis,t,u) = 0{s,t,u,p\p\p') 



(5) 



o{s',t',u',p'ip'lp'i)^d{s\t'y) 



\0{s,t,u)\ 

\ \d{s',t',u')\ 



(6) 





(YC^Y^-TE-TP)[Os+p] 



(7) 



where the function O is defined by the relation O = \0\0. 
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We can provide more motivation for tlie prescription in (6) as follows. If we 
were to drop the four terms depending on Ot in the equation for Os+p, we would 
have a much simpler equation in which it is possible to set s = and in which the 
solution depends ont^u only. In fact, if we make the replacements Os+p S and 
(t + u) then we would be left with the linearized SD equation for the two-point 
function used to study the dynamical mass S(p^). The critical coupling would be 
the same, since we have the two diagrams A and B in Fig. 5 which give identical 
one-gauge-boson exchange contributions in this case, and a factor of 1/2. 

In our case of constant coupling this equation is well-known to possess a solution 
of the form T.ijP') k, 1/|9 up to logs. We therefore deduce that the result for Og+p = 
\0s+p\0s+p is such that it falls like l/{t -\- u). This agrees with the fact that we 
are here summing the set of ladder graphs generated by diagrams A and B, and 
therefore the solution should behave like iKjP'Y k, Ijp^ when large momentum flows 
though both halves of the diagram and s = 0. We see that this consistency is a direct 
result of the prescription in (6), and we therefore proceed to use this prescription in 
the full set of equations. Note that the above discussion provides a naive guess as 
to the size of the critical coupling; it will be of interest to see how this survives in 
the full set of equations. 

5 Numerical Analysis 

We discretize the three-dimensional argument space spanned by s, t, and m, and 
the four-dimensional integration momentum space. We take the integration 4- 
momentum to have the form = {kQ,\k\ sin 6 cos ^, sin (9 sin ^, |A;| cos 6). Each of 
the variables s, t, u, |A;p and are spaced logarithmically according to log^o(^//i) + 
^logio(^[/y/^/fi)) where i = 1, ...,n. We discretize the integration angles by the pre- 
scriptions = 2Tri/n and cos^(z) = — 1 -|- 2i/(n + 1). We do not take the points 
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^ = 0, TT as integration points, because the kernel lias an integrable (coUinear) singu- 
larity there. The ratio of the ultraviolet and infrared cutoffs, Auv/^ir, should not 
be too small, but we cannot take it too large either because of the limited number 
of points on the lattice. 

Since we have specialized to pj = for alH = 1, 4, we may take a reference 
frame in which 

Pi = (Po,-|p|,0,0) 

P2 = (po,|p|,0,0) 

Pa = (po, IpI cos6', IpI sin6l, 0) 

P4 = (po, cos6l, sin6',0) (8) 

In this frame, s = ipl, u = 4|j9p(l + cos6') and t = 4|pp(l — cos^). For each 
point on the s, t, u lattice we calculate po, \p\, and 0, and then use these values 
to determine the kernels of the integral equations and the A;-dependent arguments 
s\t\u^ (i = A, ...,F) of the 4-point functions. Note that the 4-point functions are 
only defined for certain values of their arguments s, t, u according to the initial 
discretization. However, these 4-point functions are needed at different values of 
their arguments inside the integrals, due to the loop momentum. The values of the 
4-point functions at these points are found by linear interpolation. 

The relaxation method used consists of inserting initial configurations for the 
4-point functions Os+p and Ot, and then iterating (4) until it is satisfied to a 
reasonable accuracy. For the results we present, the number of points on which the 
4-point functions are defined is 12^ and the number of integration points is 12^. 
Moreover, we choose A[/y/A/^ = 6. The general features of the 4-point functions 
do not vary in any drastic way either by increasing the total number of points from 
6^ to 12^, or by varying the Af/y/A/^ ratio. 
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6 The results 



Our equations exhibit critical behavior in the sense that the coupling must be suffi- 
ciently large to support a nonvanishing solution. For Auv/^ir = 6 we find a critical 
coupling 

^^0.96 ±0.1. (9) 

etc 

We denote by a^^* the critical coupling required for a chiral symmetry breaking 
solution in the 2-point function. To be consistent, we calculate a^P* using our 4- 
point function algorithm by keeping only diagrams A and B, setting s = 0, and 
using the same A^/y/A/^. By our previous discussion this is equivalent (up to our 
approximations) to the 2-point SD equation. The ratio is not considerably 
affected by varying the kuv/ ^ir ratio. Our estimated error in etc refiects only the 
effect of the finite number of points used, which limits the accuracy of the integration 
and the interpolation procedure used inside the integrals. Rather surprisingly, it 
appears that our naive guess above as to the size of the critical coupling survives in 
the full equations. 

The numerical solution of (4) is shown in Figs. 7-12. In Figs. 7-10 we give the 
functional dependence of the function Og+p = \0s+p\0s+p on t = t/^/s on 
u = u/Ajjf^, for two different values of s = s/A|^. We use a logarithmic scale and 
only show points which are, in absolute value, greater than 10^''. Note that for 
u < t, which in our frame of reference corresponds to angles 9 > 7r/2, the 4-point 
function is negative. In that region we plot the function —Os+p, as shown in Figs. 8 
and 10. In Figs. 11-12 we give the analogous results for the function Ot = \0t\0t- 
This function is always positive, although it comes close to vanishing for large t and 
u. 



*Note that the perspective on these figures is rotated by tt with respect to Figs. 7 and 9. 
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One of the main features of the 4-point functions is their slow fah with s. 
Another is that Os+p is roughly antisymmetric and Or roughly symmetric with 
respect to the interchange t u. These symmetry properties are quite consistent 
with the form of the full integral equations in (7). When the solutions are inserted 
into these equations we see that all terms in the Os+p equation are antisymmetric 
under the interchange t u, while all terms in the Ot equation are symmetric. This 
follows from the way the C and D diagrams are related to the E and F diagrams 
when t u. Note that a solution with Os+p symmetric and Ot antisymmetric 
under t u would also be consistent, but it is disfavored by the equations. 

The order of magnitude of |Os+p| and \0t\ is the same when the angle 6 is not 
too close to 7r/2. This means that the ratios of the 4-point functions multiplying the 
coupling terms in (7) do not have a major effect; note that the appearance of these 
ratios is the main change from the equations in (4). In certain regions in the plots 
the 4-point functions display an almost linear slope (in logarithmic coordinates). 
Both IO^+pI and Or fall approximately as while remaining almost constant 
with u Avhen t < u, and fall approximately as l/u"^ while remaining almost constant 
with t when t > u. The plots also display some structure (ridges etc.) on a smaller 
scale. We believe that these latter features are artifacts of our various truncations 
and/or the discretization. 

We would like some idea as to the effect of our procedure of dropping terms 
proportional to the external momenta. One check is obtained by solving the system 
of integral equations when the terms coupling the two equations with each other 
are dropped. We plot the corresponding result for Os+p in Fig. 13, and we find a 
result which is falling roughly as l/(t u) except for the highest values oi t u. 
This behavior over most t + m is consistent with our discussion above showing the 
correspondence of this case to the two-point function. The omission of terms in the 
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integral proportional to the external momentum causes a distortion in the solution 
only at the high end oit + u. 

By comparing Fig. 13 with our full result for Os+p, we can see that the terms 
which couple the two full equations together have a very large effect on the form of 
the solution. Most dramatic is the very strong dependence of Os+p ont — u. In the 
case of Ot there is less difference between full result and the decoupled equation 
result displayed in Fig. 14. Also note that in the decoupled equation for Os+p, the 
variable s enters in a simpler way and we are able to set it equal to a constant (the 
infrared cut-off in Fig. 13). This is not possible in the full equations. Nevertheless, 
we find a fairly weak dependence on s in the full results. 

We have also made the following tests. When external momentum terms are 
dropped, we are distorting the integrand for loop momenta of order than or less 
than the external momenta. We may thus try making other distortions of the full 
equations in the infrared loop-momentum region and see the consequences. As one 
distortion, we apply an infrared cutoff on the loop momentum at a k'^ equal to the 
in (5). As another distortion, Ave replace all appearances of {P — kY, where P is one 
or some combination of the external momenta, by the quantity max(P^, k"^), which 
is angle independent. In both cases, we find that the resulting 4-point functions are 
similar to those presented above, except for some of the highest values of s, t, and 
u. 

7 Conclusions 

In this work, we have treated a linearized and truncated Schwinger-Dyson equation 
for the 4-point function of fermions interacting via a strongly interacting U(l) gauge 
interaction in the Landau gauge. We considered chirality-changing 4-point functions 
which receive nonperturbative contributions exclusively. Our basic philosophy in 
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dealing with the full intractable problem is to make various truncations to it which 
allow solutions to be found, and then identify those features of the solutions which 
are generic. This work represents a first step in this direction. 

Thus far, our numerical results indicate that the equations exhibit a critical 
behavior, and that the critical coupling appears to be roughly equal to the one 
required for the formation of 2-fermion condensates. Our results also provide a first 
indication of the momentum dependence of the 4-point functions. Further studies 
and possible ramifications will be discussed elsewhere. 
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Figure 1: The schematic form of the SD equation for the 2-point function. 



Figure 2: The schematic form of the SD equation in ladder approximation. 
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Figure 4: a) A term in the effective potential, b) An example of a leading term at 
first order in a, when the fermions are massless. c) The implied approximation. 
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Figure 5: The schematic form of the SD equation. We have labeled the four fermions 
by their spinor indices. We label the diagrams in the right-hand side by the capital 
letters A, ...,F. We omit the factor of | multiplying the right-hand side. 
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Figure 6: a) The first of the vacuum terms in the effective potential which would 
lead to a non-linear SD equation, b) The corresponding contribution to the SD 
equation. 
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Figure 7: The 4-point function Os+p when Os+p > 0, for s 




Figure 9: The 4-point function Os+p when Os+p > 0, for s = 25. 




Figure 10: The 4-point function —Os+p when Os+p < 0, for s = 25. 
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Figure 13: The 4-point function Os+p for s = 1, when the terms coupling the two 
equations are dropped. 
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